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Abstract. We construct explicit representations of the Heisenberg-Weyl algebra 
[P, M] — 1 in terms of ladder operators acting in the space of Sheffer-type polynomials. 
Thus we establish a link between the monomiality principle and the umbral calculus. 
We use certain operator identities which allow one to evaluate explicitly special boson 
matrix elements between the coherent states. This yields a general demonstration of 
boson normal ordering of operator functions linear in either creation or annihilation 
operators. We indicate possible applications of these methods in other fields. 



1. Introduction 

The Heisenberg-Weyl algebra satisfying the commutation relation 

[P,M] = 1 (1) 

has attracted attention of physicists and mathematicians for a long time. In Quantum 
Mechanics the commutator of momentum and position operators provides the most 
famous example of Eq.(£Q). Also, in the second quantization method the boson creation 
a) and annihilation a operators satisfy [a,(v] = 1. 

Among the infinity of possible representations of Eq.Q the simplest one is in terms 
of X, the multiplication by x and the derivative D = 4-, as they satisfy 

[D,X} = 1. (2) 
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(3) 



Evidently X and D acting on the monomials x n give 

Xx n = x n+l^ 

Dx n = nx 

In this note we shall be concerned with constructing and studying the representations of 
Eq-O in terms of the operators M(X, D) and P(D) such that the action of M and P on 
certain polynomials s n (x) is analogous to the action of X and D on monomials. More 
specifically we shall search for M(X, D) and P(D) and their associated polynomials 
s n (x) (of degree n, n = 0, 1, 2, ...) that satisfy 

Ms n (x) = s n+ x(x), , . 

Ps n (x) = ns n _i(x). 

The polynomials s n (x) are then called quasi-monomials with respect to M and P. These 
operators can be immediately recognized as raising and lowering operators acting on 
the s n (x)'s. Obviously M and P satisfy Eq.(JTJ. Further consequence of Eq.(JU) is the 
eigenproperty of MP 

MPs n (x) =ns n (x). (5) 

The polynomials s n (x) are obtained through the action of M n on Sq(x) 

s n (x) = M n s (x) (6) 

(in the following we shall always set sq(x) = 1), and consequently the exponential 
generating function (egf) of s n (x) is 

oo \n 

G(X,x) = J2s n (x)^- = e XM l. (7) 

71=0 n - 

Also, if we write the quasimonomial s n (x) explicitly as 

n 

s n{x) = S n,kX k , (8) 
k=0 



then 

n 



(x) 



^ Sn,kX k 
k=0 



1. (9) 



Several types of such polynomial sequences were studied recently using this 
monomiality principle first formulated in [1 j, and embodied in Eqs.(@J-(J7j). Among the 
polynomials encountered in Quantum Mechanics, Hermite and Laguerre polynomials 
are of Sheffer-type, whereas Legendre, Jacobi and Gegenbauer polynomials are not. 

Here we shall show that if s n (x) are of Sheffer-type then it is always 

possible to find explicit representations of M and P. Conversely, if M = M(X, D) and 
P = P(D) then s n (x) of Eq.(HJ) are necessarily of Sheffer-type. 



2. Sheffer-type polynomials and monomiality 



The properties of Sheffer-type polynomials are naturally handled within the so called 
umbral calculus |3J,[I]. Let us recall some relevant facts about them with special 
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emphasis on their ladder structure. Suppose we have a polynomial sequence s n (x), 
n = 0, 1,2, ... (s n (x) being a polynomial of degree n). It is called of a Sheffer A-type 
zero |2| (which we shall call here Sheffer-type) if there exists a function f(x) such that 

f{D)s n (x) = ns n -i(x), (10) 

which is the lowering operator. This characterisation is not unique, i.e. there are a lot 
Sheffer-type sequences s n (x) satisfying Eq. (fTU|) for a given f(x). We can further classify 
them by postulating the existence of the associated raising operator. A general theorem 
|3j,|lj states that a polynomial sequence s n (x) satisfying the monomiality principle 
Eq.(£H) with an operator P given as a function of the derivative operator only, is uniquely 
determined by two functions f(x) and g(x) such that /(0) = 0, /'(0) ^ and g(0) ^ 0. 
The egf of s n (x) is then equal to 

oo \n i 

V s (x)— = - e x/_1(A) (11) 

h n{) n\ g{f- l {\)) ' 

and their associated raising and lowering operators of Eq.(JU) are given by jlj 

P = f(D), 
M 



y 9'{D) 



_j_ (12) 

f'(D) ■ 

Observe the important fact that X enters M only linearly. Also, the order of X and D in 
M matters. The above holds true also for f(x) and g(x) which are formal power series. 
Any pair M, P from Eq. (jl2j) automatically satisfies Eq.Q. Here are some examples of 
so obtained representations of the Heisenberg-Weyl algebra: 

a) M(X, D) = IX — D, P{D) = \D, 

Sn{x) — H n (x) - Hermite polynomials; G(X,x) = e 2Xx ~ x2 . 

b) M(X,D) = -XD 2 + (2X - 1)D - X - 1, P(D) = - En=i D ^ 

Sn( x ) = n\L n (x) - where L n (x) are Laguerre polynomials; G(X,x) = yta^ 1 " 1 • 

c) M(X,D)=X T ± B , P(D) = -±D 2 + D, 



s n (x) = P n (x) - Bessel polynomials (THj; G(X,x) = e x(1 y/1 2X] . 

d) M{X,D)=X(1 + D), P(D)=ln(l + D), 

s n ( x ) — B n (x) - the exponential (Bell) polynomials; G(X,x) = e 3 ^ -1 ). 

e) M(X, D) = Xe~ D , P{D) = e D - 1, 



r (a+i"-n) " ^ ne l° wer factorial polynomials [TT| : G(A,x) = e xln( - 1+x \ 



f) M{X, D) = (X - tan(.D)) cos 2 (L>), P(D) = tan(D), 

s n (x) = R n (x) - Halm polynomials [12J; G(X,x) = e xarctan ( A ). 

g) M(X,D)=X 1 -^ffiD, P(D)=W L (D), 
where Wl(x) is the Lambert W function |13j : 

s n(%) = In{ x ) ~ f ne idempotent polynomials [T3j; G(X,x) = e xXe * . 
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h) The monomiality principle has also been applied to the polynomial sequences 
of more than one variable PQ,[5]-[II]. A case in point are Hermite - Kampe de Feriet 
polynomials of index m = 1, 2, of two variables x, y, defined as: 

[n/m] n-mr.r 

H^\x,y) =n\ y (13) 

n y ,yj (n-mr)\r\ y ' 

for which (here D x = ^) 

P = D X , M = x + myD™- 1 . (14) 

Their egf is a two-variable generalization of G(X, x) of example a) and it reads 

oo \ n 

G(X,x,y) = J2Ht\x,y)- = e^ A "\ (15) 
For other examples see e.g. jB]. 

3. Boson normal order and monomiality 

Now we give an application of the above formalism exploiting the analogy with the 
second quantization approach. We shall consider boson creation a) and annihilation a 
operators satisfying [a, a)} = 1. There are two sets of states of great importance in that 
representation. First, there are number states \n), (n = 0, 1, 2, ...), which are eigenstates 
of the number operator, a^a\n) = n\n), (n\n ) = 5 nn '- The action of boson operators in 
the number states basis is a^\n) = \/n + l\n + 1) and a\n) = \/n\n — 1). The second set 
are coherent states \z) = e~'^' 2//2 J2^Lo l n )' ^ ne eigenstates of the annihilation operator 
a, a\z) = z\z) [T^j. Here we will show how to use the above formulas to calculate some 
special matrix elements of operator functions, which are linear in either o or at As a 
byproduct we shall produce relations for the normally ordered form of these operator 
functions. 

To make the analogy with the monomiality property Eqs.(jnj) and (j3J), it is convenient 
to redefine the number states 



\n) = Vn\\n) (16) 

(note that |0) = |0)). Then the creation and annihilation operators act as 

a t|n) = \n+l) (17) 
a\n) = n\n — 1) (18) 

which allows us now to make the correspondence: 
X < ► at 

D < ► a (19) 

x n < — > \n) , n = 0,l,2,... . 

Then the operators M and P take the form (see Eq.([12|)) 

P(a) = f(a), 



M(a, a 



at - 2^ 

9{a) 



1 



(20) 



/'(a) • 
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Recalling Eqs.®,© and © we get 

n 

M(a,at)]>) = ]>> n , fc (at) fc |0>. (21) 

k=0 

In the coherent states representation it yields 

(z\ [Af (a, a f )] n |0) = s n (z*)(z\0). (22) 
Exponentiating M(a, a)) and using Eq. (jll|) we obtain 



(z\e 



AAf(o,ot)|Q\ 



z\0). 



(23) 



By the same token one obtains closed form expressions for the following matrix elements 
( |Z) is the Z-th number state, I = 0, 1, 2, ...): 



r x i n \ 

(z\ [M{a,J)\ |Z) = — Sn+l (z*)(z\0). 



and 



(z\e XM ^ at) \l) 



1 d l 

~7T\d\ l 



(z\0). 



(24) 



(25) 



[g(f- l W) 

The result of Eq. (j2*3*|) can be further extended to a general matrix element (z\e XM( - a,a ^\z') . 
To this end recall ^3] that \z) = e -l z l 2 / 2 e zat |Q) and write 



(z\e 



\z) = e 

= e-\^l'\z\e- < 

_ z*z'-\z'\ 2 /2 



10} 



-\z'\ 2 /2/ v \ e z'a! p -z'at e \M(a,rf) e z'cJ 
z'at e \M(a,ai) e z'ahQ\ 



(z\e 



(z\e 

'at = 

WHa,r)in = e z*z'~\z'\ 2 /2^ e \M(a+z'^)^ 



Next, using the property e~ z ' a M(a, a j )e z ' a = M(a + z', a)) [To] we arrive at 



> ( + g'(a + z') \ 1 

2 '~ |z ' |2/2 (z|e V ~ 9(a+z,) / f ' (a+z,) 10). 



Now we are almost ready to apply Eq.(|23|) to evaluate the matrix element on the r.h.s. 
of the above equation. Before doing so we have to appropriately redefine the functions 
f(x) - />) = f(x + z>)-f(z') and ^(s) - = g{x + z!)lg(2>). Then /(0) = and 
g(0) = 1 as required by Sheffer property for / and g. (Note that these conditions are 
not fulfilled by f(x + z') and g(x + z').) We can now write 



\(n* B'(*+z') \ 1 , / t §'(») 

/ z |g V 9 (a + Z ') J /'(a + z') |Q^ _ V «(") / . 



10) 



By going back to the initial functions / and g this readily gives the final result 



(z\e 



XM(a,at) I J 



g{f-\\ + f{z'))) 



^z*ir\x+f(z'))-z'] 



(26) 



* / 1 I ' 1 2 1 I 1 2 

where (zjz ) = e z z ~2< z I ~2l 2 l is the coherent states overlapping factor. 
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It is natural to make contact now with our solution of the normal ordering problem 
[T7j . For a general function F(a, at) its normally ordered form M F(a, a)) = F(a, at) is 
obtained by moving all the annihilation operators a to the right, using the commutation 
relations. We may additionally define the operation : G(a, at) : which means normally 
order G(a, at) without taking into account the commutation relations. Using the latter 
operation the normal ordering problem is solved for F(a, at) if we are able to find an 
operator G(a, at) for which F(a, at) =:G(a, at) : is satisfied. 

To obtain normally ordered form of e AM(a,at ) we apply the crucial property of the 
coherent state representation (see ^3] and JEj). This is that if for an arbitrary operator 
F(a, at) we have 

(z\F(a,at)\z') = (z\z') G(z*,z') (27) 
then the normally ordered form of F(a, at) is given by 

tf[F(a,at)] = : G(at,a) : . (28) 
Eqs. fl26|) and ()28|) then provide the central result 



3 AM(a,at) 



,at[/-l(A+/(o))-a]. 



9W 



(29) 



g(f-i(\ + f(a))) 

Let us point out again that at appears linearly in M(a,at), see Eq.(pUj). By hermitian 
conjugation of Eq. (|29|) we obtain the expression for the normal form of e AMt ( a > at ) ; where 
M^(a,a^) is linear in a. In this context we refer to our previous work where 
the results equivalent to Eqs.(|2Hjl and (j2H| were obtained using the substitution group 
approach 0. 

We shall conclude by enumerating some examples of evaluation of coherent state matrix 
elements of Eqs.()22|) and (|2*S|). We choose the M(a, a^)'s as in the list a) - g) above: 

a) (z\{-a + 2at) n \0) = H n (z*){z\0), 

b) (z\ \-a ] a + (2at - l)a - at + ll" |0) 



c) 



2 | e A[-ota+(2ot-l)a-at+l] in 



1 



n\L n „i(z*)(z\0), 

Z*\- 



l-A(z'-l)' 



A(l-z')-l 



(z\z>). 



z\ a 



1-a 



\0) = P n (z*)(z\0), 
z \e x ( af Th)\ z ') = e ^-V^i^ Z W)- z ']( z \z'). 
d) (z\(ata + at) n \0) = B n (^)(z\0), 

z \ e \(a^a+a*)\ z /\ e z*(z'+l)(e A -l) / z \z'). 

,|(ate-T|0) = r gQ^(^|0), 
f) (z\ [(at -tan(a))cos 2 (a)l"|0) = R n (z*){z\0), 



j^ z | gA(at — tan(a)) cos 2 (a) |^/^ 



cos [arctan( A+tan (z') )] 
cos(z') 



,z*[arctan(Atan(z'))— z'] /^|^'\ 
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g) (z\[a^-^a] n \0) = I n (z*)(z\0), 

4. Conclusions and outlook 

In this work we have used the prescription from the umbral calculus to find explicit forms 
of the raising and lowering operators acting in the space of Sheffer-type polynomials. 
Their specific form given by Eq. (jl2j) permitted us to calculate certain coherent state 
matrix elements. Those in turn have given directly, see Eq.()29|). an explicit expression 
for the normally ordered form of operator functions of M(a, a)) since M is linear in either 
o or at A legitimate question is if other types of sequences would be quasimonomials 
with respect to some operators M and P. The affirmative answer was given in [TH] 
where it has been demonstrated that for any sequence such operators can be found. 
For a number of polynomial sequences the operators M and P can actually be written 
down explicitly. However, only in the Sheffer case P = P(D) and M(X, D) is linear in 
X. In any other case P = P(X, D) and and it appears that this circumstance renders 
the calculation of coherent state matrix elements rather difficult. This is without doubt 
another demonstration of an intrinsic difficulty to perform a normal ordering of arbitrary 
operator function. 

Let us mention that the operator methods based on monomiality principle 
developed above have interesting, and as yet non-explored, ramifications to at least 
two other fields [TH]: non-local evolution equations and generalised heat equations. We 
shall briefly sketch these points here, leaving the details for the forthcoming publication. 
Let us consider here the case c) above for which s n (x) are the Bessel polynomials with 
M(X, D) = Xj 1 ^ and P(D) = -\D 2 + D. Then it is possible to look for the solution 
of the following operator evolution equation: 

|-F(X, y) = [P(D) + M(X, D)]F(X, y), (30) 

with the 'initial' condition F(X, 0) = q(X). Using A' 1 = J °° e~ sA ds we find 
M = X J °° e- s e sD ds, and then with Mq(X) = X / °° e~ s q(X + s)ds, Eq.® takes 
the following non-local form 

(9 1 r°° 

—F(X,y) = [--D 2 + D}F(X,y)+X e~ s F(X + s,y)ds, (31) 
oy I Jo 

whose formal solution can be written as 

F(X, y) = e y{p+M) q(X) = e~^ 2 e yP e yM q{X). (32) 

For given q{X), e yM q(X) can be evaluated recursively through e yM q(X) = 
E^o^WS with 7r «( X ) = A/ °°e- s 7r n _i(X + s)ds. With additional standard 
transformations a closed form of F(X, y) can be obtained. 

Another perspective will be offered by further two variable [20] extension of Eq. ()14|) . 
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Define 

n = f(D x ), (33) 
= M x + 2M y f(D y ) (34) 

with D v = |j and M,, = — ^Jjt^ for = z,y. Evidently [11,9] = 1. This 
particular two- variable generalisation of monomiality operators allows one the following 
extension of Hermite-Kampe de Feriet polynomials for m = 2 of Eq. (jl3|) : 

[—] 

S n {x,y)=n\2^— f v 35 

f^ {n — 2r)\r\ 

which are umbral extensions |2Hj of the ordinary case, where s n (rj) are Sheffer 
polynomials associated with / and g. It is worth noting that S n (x,y) are solutions 
of generalized heat equation 

f(D y )S n (x,y) = [f(D x )) 2 S n (x,y). (36) 

These and related topics will be developed in subsequent publications. 

We thank L. Haddad for important discussions. 

PB wishes to thank the Polish Ministry of Scientific Research and Information 
Technology for support under Grant no: 1P03B 051 26. 
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